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Abstract 

A Green's function technique for a modified Ising model in a transverse field is applied, which 
allows to calculate the damping of the elementary excitations and the phase transition temperature 



\ of ferroelectric thin films with structural defects. Based on an analytical expression for the damp- 

O ' ing function, we analyze its dependence on temperature, film thickness and interaction strength 

o : 

numerically. The results demonstrate that defect layers in ferroelectric thin films, layers with im- 

^ , purities or vacancies as well as layers with dislocations are able to induce a strong increase of the 

00 ' damping due to different exchange interactions within the defect layers. The results are in good 

(N ■ 

OO ■ agreement with experimental data for thin ferroelectric films with different thickness. 
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I. INTRODUCTION 



Defects in crystals influence the physical properties of almost all materials significantly. 
In the last decades, defect engineering has been developed as a part of modern semicon- 
ductor technology. There is also an increasing interest in studying defects and their related 
strain fields in ferroelectrics, in particular in low dimensional systems Q]. The influence of 
defects on the phase transition is one of key-like topics in the context of recent studies of 
ferroelectrics. The permanent attention attracted to ferroelectric (FE) material, is related 
to its role played in real ferroelectrics concerning fundamental problems as well as the broad 
variety of applications based on thin FE films, especially due to their suitability for appli- 
cation in nonvolatile FE random access memories P]. It is well known that localized spin 
excitations can arise in FE and magnetic systems with broken spatial translational symme- 
try, i.e. due to the presence of boundary surfaces, interfaces, impurities and other defects 
in the crystal. These localized modes appearing in the above mentioned structures are ex- 
perimentally observed using the neutron and Raman scattering 0, |^ or, more recently, far 
infrared absorption jsj. From more theoretical point of view the Green's function theory has 
)een extensively employed to describe these modes in semi-infinite transverse Ising models 
611 and FE thin films Ml- In addition to bulk spin waves there may occur localized modes 

n 

associated with the impurity layer |63 and with the surface J'il]. More phenomenologically 
the influence of molecular impurity ions on FE phase transitions is studied by Vikhin and 
Maksimova using the Landau theory Qj. 

The relaxation (damping) of polarization motion is a result of microscopic fluctuations 
of the spins by means a coupling to other degrees of freedom such as phonons, impurities, 
defects etc. The mechanism behind can be originated by different reasons. Once there 
could be realized a direct coupling of the elementary spin wave excitations to a thermal 
bath or to phonons manifested by a spin-phonon scattering. Otherwise more complicated 
microscopic processes, such as "slow relaxing" impurities could control the fluctuations and 
the relaxation processes of the system. In that case the impurity modes are damped which 
leads to a relaxation of the polarization motion due to a direct coupling. In the present paper 
we want to study such a mechanism in detail based on an Ising model in a transverse field. 
There are two kind of relaxation processes: inherent (pure bulk material) and non-inherent 
(due to defects, impurities, surface, finite size). Each microscopic relaxation mechanism 
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predicts its own temperature and frequency dependencies. The damping can also depend 
on defect /impurity concentration and the sample size. Using X-ray diffraction it is shown 
that too many structural defects can produce an excessive line broadening I9I, e.g. defects 
contribute to the line width and consequently the damping increases. Doping of thin films 
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|. Ingvarsson et 



with different ions yields to variable increase of the damping 
al. showed that the relaxation mechanism in thin films could be realized in analogy to 
the bulk relaxation, where the phonon scattering process in the bulk should be replaced by 
surface and defect scattering in thin films. 

A first attempt to find out the influence of defects on the polarization of FE thin films are 



elucidated by Alpay et al. 



15|. Based on the Isin 



model in a transverse field and using 



the Green's function theory Wesselinowa et al. [ig] had been successfully in calculating 
the spin-wave energies, the polarization and the phase transition temperature for a FE thin 
film with different structural defect layers. The aim of the present paper is to study the 

damping effects and the relaxation times in FE thin films with different structural defect 

r— 1 

layers. To that purpose a previous paper [IT] is extended and modified in order to obtain 
the temperature dependence of the soft modes and their damping within FE thin films. It 
was obtained that with increasing film thickness the frequency increases likewise and the 
damping is larger for thin films than that for the bulk material. This is in agreement with 
different experimental observations Using Raman spectroscopy the temperature 

dependence of the phonon modes for thin FE films of PbTiOa is discussed by Taguchi et 
al. and Fu et al. It was shown that in comparison with the single crystal spectra, 
the Raman frequencies for the thin film are shifted remarkably to low frequencies and that 
the Raman lines are broadened. Similar results were obtained for PbTiOs fine particles by 
Ishikawa et al. 



II. THE MODEL AND THE MATRIX GREEN'S FUNCTION 

Let us consider a three dimensional ferroelectric system on a simple cubic lattice composed 
of layers in z-direction. The layers are numbered by n = 1, ...A^, where the layers n = 1 
and n = N represent the two surfaces of the system. The bulk is established by the remaining 
(A^ — 2) layers. The specific surface effects are included by additional coupling parameters 
between bulk and surface layers. In particular, we start with the Hamiltonian of the Ising 
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model in a transverse field which includes both, bulk and surface properties: 



(1) 



ij ieb ies 

where and are components of spin-| operators, Qf, and Qg represent transverse fields 
in the bulk and surface layers, and the sums are over the internal and surface lattice points, 
respectively. Jij is an exchange interaction between spins at nearest-neighbor sites i and j, 
and Jij = Js between spins on the surface layer, otherwise it is Jb- We assume that one or 
more of the layers can be defect, since Jd and Qd denote the exchange interaction and the 
transverse field of the defect layer. The ordered phase is characterized by the non zero mean 
values (S^) ^ and (S^) ^ 0. Hence it is appropriate to choose a new coordinate system by 
rotating the original one, used in Eq. (P), by the angle 9 in the xz plane [2^. The rotation 
angle 9 is determined by the requirement {S^') = in the new coordinate system. 
The retarded Green's function to be calculated is defined as 



=«5+(t);57(0)», 



(2) 



where and 5" are the spin-| operators in the rotated system. On introducing the two- 
dimensional Fourier transform Gmn O^Wj^^), one has the following form: 



5 Sj 



a 



;^exp(zk|j(ri - rj))Gn,n,m,uj), 



(3) 



where N' is the number of sites in any of the lattice planes, rj and Ui represent the position 
vectors of site i and the layer index, respectively, ky = {kx,ky) is a two-dimensional wave 
vector parallel to the surface. The summation is taken over the Brillouin zone. 

As a result the equation of motion for the Green's function Eq. (jS)) of the ferroelectric 
thin film for T < T^. has the following matrix form: 



H(a;)G(k||,a;) = R, 



(4) 



where H can be expressed as: 
/ u — Vi + iji 



H 



k2 




00 -V2 + 272 
ks uj 





k2 

^3 + ^73 k3 



/cat — V^r + i^N J 
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with 



kn = JbUnSin^On, n^l,...,N, 

Vi = sin 01 + -aiJs cos^ 9i - sin^ ^i7(k||) + JftfTa cos^ 62, 

= 20ft sin 6^2 + ^0-2 Jft cos^ 62 - sin^ ^27(k||) + JgCri cos^ 61 + Jftcra cos^ 6'3, 

K = 20„ sin 0n + ^crnJb cos^ 6'„ - sin^ 6'n7(k|| ) + J„-icr„_i cos^ 6'„_i 

Viv = 2fisSin6'jv + ^o-jv-ZsCOS^^jv — — y^sin^6'iv7(k||) + J6(77v-iCOS^6'iv-i, 
2N^ 



7n = TTicH J] [(^n(qihk|i-qii) + ^n(k|i-Pii-qihPii + qii))^ 



P||,qil 

* [nn{P\\)i(^n + finiPW + ^\\) + ^n(k|| - q||)) 

- n„(pi| + q||)n„(k|| - qy)] 

* ^(en(k|| - qy) + e„(p|| + qy) - €n(p||) - e„(k||)) 

+ [( J„_i7(q||) cos^ 9n^if + ( J„-i7(k|| - P|| - q\\) cos^ O^^if] 

* [n„_i(p||)(cr„_i + n„_i(p|| + qy) + nn-i(k|| - qy)) 

- nn_i(p|| + q||)nn_i(k|| - qy)] 

* 5(e„_i(ky - qy) + e„_i(py + qy) - e„_i(py) - e„_i(ky)) 

+ [(J„+i7(qy) cos^6'n+i)^ + (J„+i7(ky - Py - qy) cos^ 6'„+i)^] 

* [^n+i(P||)((^n+i + ^n+i(P|| + Qn) + ?in+i(ky - qy)) 

- ^n+i(P|| + qy)n„+i(ky - qy)] 

* 5(en+i(ky -qy) + e„+i(py +qy) -en+i(py) -en+i(ky)) , 
t^(q|hky -qy) = J(qy) cos^ ^ - ^ J(ky - qy) sin^ ^, 

7(ky) = i(cOs(M+COs(M)- 

Here we have introduced the notations Ji = Jn = Js, Jn = Jb for n = 2, 3, 4, iV — 1, 
Qi = Qn = Qs, 0„ = Qi)n = 2,3,4, N — 1), Jq = Jn+i = 0. The quantity cr(T) is the 
relative polarization in the direction of the mean field and is equal to 2{S^'). For the rotation 
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angle 9 we have the following two solutions in the generalized Hartree-Fock approximation: 

1. cos^ = 0, i.e. = ^, if T>Tc; 

2. sme= — = —, a T<T,. 

a J o 

In order to obtain the solutions of the matrix Eq. we introduce the two-dimensional 
column matrices, Gm and Rm, where the elements are given by (G„)m = Gmn and (Rn)m = 
o"n^mn, SO that Eq. (jH) yields 

H(a;)G„ = R,,. (5) 

From Eq. (jSJ, G'„„(a;) is obtained as: 

OU^ - (6) 

The quantity is the determinant made by replacing the ra-th column of the deter- 

minant |-f/^(a')| by Rn- The poles a;„ of the Green's function Gn„(a') can be calculated by 
solving \H{{jS) \ = 0. 

The relative polarization of the ?7,-th layer is given by 

,2 , 



a-n 



/cr„J„ ^ 1 - 0.5 sin 6'„7(k|j) ^, a;„\-i ,„■, 



Eq. ((Zj) has to be numerically calculated. Due to the assumption of symmetrical surfaces, 
there are layer polarizations, which have to be solved self-consistently. In order to obtain 
the dependence of the Curie temperature on the film thickness A^, we let all cr's be zero 
in Eq. ((7j) and solve the expressions self-consistently. 



III. NUMERICAL RESULTS AND DISCUSSION 



In this section we shall present the numerical calculations of our theoretical results taking 
the following model parameters: Jf, = 495 K, fi^ = 20 K. We have calculated the temperature 
dependence of the damping from 7 = ;^ J2n 7" ^ simple cubic thin film (ky = 0) and for 
different values of the exchange interaction constants. The numerical results expose some 
interesting and novel characteristics in the damping values in comparison to the case of FE 
thin films without defects. The results for film thickness N = 7 and different exchange 
interaction parameters in the defect layer are presented in Fig. We consider firstly 
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the case where the middle layer is defect, which is possibly the case when the layer has 
vacancies or impurities with smaller radius and larger distances between them. The exchange 
interaction is smaller than the value of the bulk interaction Jf, = 495K and has the value 
Ja = Jd = 300 K (Fig. m curve 1). The damping is larger than in the case without defects, 
J4 = Jh (Fig. m curve 3). This is in agreement with the experimental data 



Unfortunately, there are more experimental data for damping of magnetic thin films and not 
for ferroelectric thin films. The damping increases with increasing temperature, T Tc of 
the thin film. This could explain the observed experimental results of Raman scattering from 
ferroelectric thin films that the line shapes of the film become broad as the temperatures 
approaches The critical temperature decreases due to the smaller value. 

From Eq. ((Tj) by solving (j{Tc) = we get the Curie temperatures: Tc(bulk) = 219.5 K and 



Tc{Jd = 300) K = 213.8 K. Pontes et al. [2l| carried out dielectric and Raman spectroscopy 



studies and obtained that with addition of Sr to PbTiOs, the phase transition temperature 
decreases with the increasing Sr concentration. 

For the case where J4 = = 1000 K (Fig.HJ curve 2), i.e. Jd is larger than the value of the 
bulk interaction constant Jb, (for example when the impurities have a larger radius compared 
with the constituent ions) the damping is again larger in comparison to the case without 
defects, Jd = Jh- Experimentally one observes an increase of the damping in presence of 



defects ll 
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12l Il3| . The value of of the film is enhanced in comparison to the bulk 
value without defects due to larger Jd values. It results: Tc{Jd = 1000) K = 228.5 K. This 
is the opposite behavior compared to the case of Jd = 300 K, Jd < Jb- The second case, 
where Jd > Jb could explain the experimentally obtained increase in the phase transition 
temperature by the substitution of impurities, such as Nd, V, W and Nb in the layers of 



BTO films 2^ 



3, 



2J| or by increasing Ca contents in SET thin films j2 



We obtain that the critical temperatures of the FE phase transition are increased or 
decreased due to different exchange interactions in the defect layers. Our results are in qual- 
itatively agreement with the experimental data of Noguchi et al. [2^. They have studied 
SET films {Tc = 295° C) by the substitution of rare earth cations of La, Ce, Pr, Nd and 
Sm as well as Ei at the A site (Sr site) with Sr vacancies and have shown, that La modifi- 
cation induces soft behavior (lower and lower Pr), while a large amount of Nd and Sm 
substitution results in a very high (hard), as a result of defect engineering of both Sr 
and oxide vacancies. For not only SET, but also other bismuth layer-structured ferroelec- 
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trie (BLSF) Tc is strongly influenced by of A-site cations, and BLSF with smaller A-site 
cations (Ca^^) tend to show a higher Tc (420° C). The same amount of larger Ba^+ brings 
about a relaxation of FE distortions and leads to a decrease of to 120*^C. The substitution 
of La led to a marked decrease in to 180° C (a;=0.5), because the induced A-site vacancies 



weaken the coupling between neighboring BOg octahedra j26|,|27|. This result corresponds 



in our calculations to the case of smaller values of the interaction constant in the defect layer 
Jd < Jb- For La-modified PbTiOs, Tc decreased significantly, too, with an increase in La 

n 

content |23]. For Bi-SBT Tc rose strongly to 405 C (x=0.2) 25]. The increase in Tc by Bi 
substitution is an opposite tendency to that in the case for La-SBT. The dominant role plays 
here the bonding characteristics with oxide ions. The influence of the orbital hybridization 
on Tc is very large, and Bi substitution resulted in a higher Tc 2^. This experimental result 
can be described qualitatively good using in our model exchange interaction parameters for 
the defect layer > Jb- 

The damping depends also on the number of inner defect layers. This is shown in Fig. |2l 
for a thin film with N = 7 layers. It is demonstrated that for Js = 300K with an increasing 
number of defect layers (curve 1 - J3 = J4 = J5 = = 600K, curve 2 - J4 = Jd = 600^) 
the damping increases, i.e. we obtain a line broadening. This is in agreement with the 
experimental data of BouUe et al. who argued that too many structural defects can 
produce excessive line broadening. 

In Fig. IHlwe have studied a ferroelectric thin film with different film thickness where one 
layer, the middle one is defect, Jd = lOOOK (A^ = 7 - curve 1) and (A = 9 - curve 2). 
The damping increases with decreasing of the film thickness, what is in agreement with our 



previous calculations without defects 



17| and with the experimental data |l8, 



IS, 2fl. The 



thinner the films the larger is the damping. 



IV. CONCLUSIONS 



Based on a modified transverse Ising model and using a Green's function technique, the 
damping and the phase transition temperature for ferroelectric thin films with structural 
defects are calculated for the first time. The dependence on temperature, film thickness and 
interaction constants is discussed. It is shown that defect layers in FE thin films, layers 
with impurities or vacancies, as well as layers with dislocations, can induce strong increase 
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or decrease of the critical temperature of FE phase transition due to different exchange 
interactions in the defect layers. We obtain that the damping in the thin films is always 
larger in comparison to the corresponding bulk values or to films without defects. There 
are different mechanisms in thin films which contribute additive to the damping, due to 
spin-spin, spin-phonon, surface, defects scattering etc. The damping in thin films is greater 



compared to the bulk case for example due to surface effects [17[ then due to spin-phonon 
interactions [29] and to defects which are considered in the present paper. The damping 
is related to the width of the half-maximum for the Raman scattering lines. It can be 



concluded, that if we want to obtain and to explain the 



experimentally in thin films and nanoparticles 



19 



arg e line broadening effects observed 
20| all contributions to the damping 



must be taken into account and that the defects play an important role. 
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T[K] 

FIG. 1: Temperature dependence of the damping 7 in cm^^ for a FE thin film with J5 = 495 K, 
= 20 K, Js = 900 K, ns = nh, N = 7 and different J^- values: (1): Ja = 300 K; (2): 1000 K; (3): 
495 K. 
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FIG. 3: Temperature dependence of the damping 7 in em ^ for J5 = 495 K, 17;, 
600K, Jd = 1000K,J7s = and different film thickness: (1): N = 7; (2): N = 9. 
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